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Abstract 

Crossover quintessence predicts that the time variation of fundamental 
couplings is substantially faster at redshift z ~ 2 than over the past few billion 
years. This could reconcile the reported time variation of the fine structure 
constant from quasar absorption lines with severe bounds from archeo-nuclear 
physics and high precision tests of the equivalence principle. We present a 
model that is consistent with all present data on time varying couplings, tests 
of the equivalence principle and cosmology. 
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Quintessence or a time varying dark energy can be associated with a scalar 
field whose cosmological expectation value varies during the recent history of the 
universe [1]. Generically, this scalar field - the cosmon - may also couple to matter 
and radiation. As a consequence, the values of "fundamental constants" like the 
fine structure constant or the ratio between the nucleon mass and the Planck mass 
depend on the value of the cosmon field and therefore on cosmological time [l]-[2]. 
Recently, a time variation of the fine structure constant has been reported [3] from 
the observation of quasar absorption lines (QSO) at redshift z 2. A typical value 
corresponds to 

A ^ = 2) =-7.10- 6 . (1) 

"em 

Similar observations infer 2 a substantially smaller value of |Aa em | at lower redshift 
z < 0.7. Furthermore, from the Oklo natural nuclear reactor one obtains a typical 
bound [5], [6] 

\Aa em (z = 0.13)| < 1Q _ 7 ^ 

" em 

and analysis of the decay rate Re 187 — > Os 187 from meteorites restricts [6] 

\Aa em (z = 0^45)1 <3 . 1Q _ 7 _ (3) 



As an obvious question one may ask if an increase of |Aa em | by almost two 
orders of magnitude between z = 0.13 and z = 2 can reasonably be explained by 
quintessence models. In this note we demonstrate that this can indeed be the case 
in a class of models of "crossover quintessence" proposed recently [7]. This contrasts 
to constant rates dAa em /dz = const, or dAa em /dt = const, for which the values 
and bounds (l)-(3) clearly are in discrepancy. The deviation from constant rates 
also affects strongly the comparison of eq. (1) with bounds from nuclear synthesis 
[8], [9] and the cosmic microwave background (CMB) [9]. 

We are aware that the bound from the Oklo natural reactor is subject to sub- 
stantial uncertainties due to possible cancellations between effects from the variation 
of a em and other "fundamental constants" like the mass ratio between pion and nu- 
cleon mass m n /m n or the weak interaction rates. Also the QSO observations need 
further verification and investigation of systematics. Nevertheless, it is well con- 
ceivable that we can get a reasonable "measurement" of the function Aa em (z) in 
a not too distant future. Furthermore, the derivative dAa em (z)/dz at z = can 
be related 3 [10] to precision tests of the equivalence principle [11]. The aim of this 
letter is to present a sample computation how the knowledge of Aa em (z) can be 
used to constrain models of quintessence. For this purpose we will use the value (1) 
as a "benchmark value" which fixes the strength of the cosmon coupling to matter 



2 A recent collection of results on Aa em as well as the time variation of other fundamental 
constants can be found in [4] . 

3 See [2] for a different relation which gives a much smaller violation of the equivalence principle. 
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and radiation. The time history Aa em (z) can then be related to the time history of 
quintessence. 

In particular, we may investigate the ratio 



A bound R < 1/50 strongly favors [10] quintessence with a time varying equation 
of state Wh = Ph/ph, where the value of (1 + Wh) at present is substantially smaller 
than for z — 2. A logarithmic dependence of Aa em on the cosmon field x leads to 
a bound for the present equation of state [10] 



In this note we investigate how the detailed dependence of a em on \ influences the 
time history of Aa em (z). We argue that reasonable /9-functions for the "running" 
of a grand unified coupling with ln% may lead to nonlinearities in a em (lnx). Those 
are quantitatively important, influencing R within a factor two. The overall picture 
remains rather solid, however: a small value of \R\ requires a quintessence model 
where the evolution of the cosmon field has considerably slowed down in the recent 
history as compared to a redshift z »s 2. This feature is characteristic for crossover 
quintessence and will not be shared by arbitrary quintessence models. In particular, 
a constant equation of state with Wh independent of z over a range < z < 3 will 
have severe difficulties to explain R < 0.1. This demonstrates how a measurement of 
Aa em (z) could become an important ingredient for the determination of the nature 
of dark energy. 

Crossover quintessence (CQ) can be characterized by a recent increase of the 
fraction in homogeneous dark energy Qh = Ph/ Per from a small but not negligible 
value for early cosmology to a value Qh ~ 0.7 today. A small early value (say 
Qh = 0.01) typically results from a cosmic attractor solution [1, 12, 13] ("tracker 
solution" ) independently of the detailed initial conditions. For such models of "early 
quintessence" the dark energy is diluted in early cosmology at the same pace as the 
dominant energy component. Dark energy and radiation or matter have therefore 
always been of a comparable magnitude and excessive "fine tuning" can be avoided 
[1, 14]. The amount of early quintessence is accessible to observation by nucleosyn- 
thesis [1, 15], CMB [16] or structure formation [17, 18]. 

The crossover to domination by dark energy must have happened after the forma- 
tion of structure. On the level of the equation of state CQ corresponds to a crossover 
from a value Wh ~ 1/3 during the radiation dominated universe to a substantially 
negative value for z < 0.5. During the crossover period the time variation of the 
cosmon field remains substantial - this may extend to z ~ 1.5 and be the origin of 
the QSO observation (1). More recently, towards the end of the crossover period, 
the equation of state may have approached the value Wh ~ — 1, corresponding to a 
suppressed cosmon kinetic energy T = (1 + Wh)ph/2 and therefore a slowing down 
of the evolution of the cosmon. For a present value 1 + <C 1 the evolution of 



R = Aa em (z = 0.13)/Aa em (z = 2). 



(4) 



wf ] < -0.9. 



(5) 
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the cosmon and the fundamental couplings has almost stopped, explaining the small 
bound (2). We note that the fraction of energy density in matter, f2 m , decreases 
rapidly once Wh has come close to —1 at the end of the crossover. The observed 
present value k, 0.3 therefore indicates that the effective "stop" in the cosmon 
evolution should have happened rather recently - a value z ~ 1 for the drop in the 
variation rate of the fundamental couplings seems rather natural in this respect. 

In this letter we want to give a specific example for this scenario. For this 
purpose we elaborate on a recent proposal [7] that CQ is connected to the existence 
of a conformal fixed point in a fundamental theory. The vicinity to a fixed point 
can give a natural explanation why the time variation of fundamental couplings is 
only a tiny effect. Our model resembles in this respect the "runaway dilaton" in 
string theories [19]. We want to demonstrate here the existence of a model which 
is compatible with (l)-(3) rather than to advocate the detailed choice of the (3- 
functions for the running couplings below. 

Let us assume that for low momenta (q 2 <C x 2 ) the yet unknown unified theory 
results in an effective grand unified model in four dimensions. There the coupling of 
the cosmon field x to gravity and the gauge fields is described by an effective action 

S = Jd 4 x^^x 2 R+l(S(x)-^xd,X 

+m\ 2 + ^F^F,A. (6) 
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We note that the value of the cosmon field x plays the role of a dynamical mass scale 
in the fundamental theory. Typically, it could be associated with the characteristic 
scale for the transition of a higher dimensional theory to an effective four dimensional 
"low energy" theory. We have normalized x sucn that its present value equals the 
reduced Planck mass M,M 2 = 1/(8tvGn). Typical particle masses are ~ x- Up to 
the x-dependence of the running dimensionless couplings 5 and Zp and the explicit 
mass term m 2 x 2 the effective action is dilatation invariant. The global dilatation 
symmetry x —> ®X is spontaneously broken by any nonzero cosmological value of X- 
The running of the effective grand unified gauge coupling g\ = g 2 /Z F is deter- 
mined by the x-dependence of Z F . The assumed fixed point behavior for a x = 
g x /4:7i is cast into the form 

da ( S \ 

P a = b 2 a x - b 4 a x - b 6 a 3 x I 6 - J . (7) 



d In x \ 1 + J a 5 

We take here a small value of |6 6 | and impose a fixed point value 

« ^ « 45 (8) 

as extrapolated from the unification of the observed gauge couplings of the standard 
model. The term ~ b 6 accounts for the dependence of the fixed point on the value 
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of the other dimensionless coupling 5. As mentioned above, we do not assign any 
particular significance to the specific form of (5 a . We only assume the existence of a 
fixed point in the approximation of constant 5 which depends slightly on the value 
of 8. In our case the approximate shift in the location of the fixed point between 
5 = and 5 = oo amounts to 4 Aax,* = b$a x J (&2^a)- Many other possible shapes 
of (3 a will lead to similar conclusions provided that one can give a reason why the 
variation of «i is small. 

For b 2 > and cxx(x) > a x,* the grand unified gauge coupling increases with 
decreasing x- Similar to QCD it grows large at some nonperturbative scale ("con- 
finement scale") that we associate with m. Our model therefore has an "infrared 
scale" m which is generated by dimensional transmutation from the running of the 
grand unified gauge coupling. We assume that the non-perturbative scale m de- 
termines the mass of the cosmon. Neglecting the effects of the explicit dilatation 
symmetry breaking reflected by m the cosmon potential would be flat. 

On the other hand the gauge coupling runs towards its ultraviolet fixed point 
(8) as x increases. At the fixed point the effective action would have an effective 
dilatation symmetry if the running of S(x) could be neglected. We will assume, 
however, that the running of 5 is unstable towards large x an d postulate (E > 0) 

m m = ^r-, (»> 



<91nx 1 + Js5 

For small enough Si = S(x = tu) the dimensionless coupling 8 stays small for a large 
range m < x < Xc- For x ^> m we may neglect the small mass m for the evolution 
of the dimensionless couplings. The evolution in this range is then governed by the 
vicinity to a conformal fixed point [7] at S — S* — , ax = «!,*■ In the same 
approximation our model corresponds to a flat direction 5 in the effective cosmon 
potential. All cosmon interactions are therefore given by its derivative coupling to 
the graviton and gauge fields. 

Nevertheless, as x increases further an effective "ultraviolet scale" Xc character- 
izes the region where 5 grows large according to the approximate solution 6 

«*> * TO)- (10) 

The ratio between the ultraviolet and infrared scales turns out exponentially huge 



4 We have introduced J a such that the shift remains finite for S — > oo. 

5 In string theories such flat directions are often called "moduli". In presence of additional 
matter fields, e. g. scalars and fermions, the flat direction in the space of scalar fields corresponds 
to constant ratios of all masses [7] , up to the small effects from varying couplings discussed in this 
note. The flat direction is consistent with the conformal symmetry, but not required by it. See [7] 
for a more detailed motivation of a flat direction in the presence of quantum effects. 

6 We take a small value J$ = 0.05 such that only the behavior for large S is affected by Jg. The 
cosmology depends only very little on J$ in this case. 
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At the crossover scale \ c the flow of the couplings witnesses a crossover from the 
vicinity of the conformal fixed point (small 5) to another regime for large 5 - in our 
case a fixed point for the effective cosmon-gravity coupling at d^ 1 = 0. We will 
see that the crossover behavior in the present epoch of the cosmological evolution 
precisely corresponds to the crossover from small to large 5 for x i n the vicinity of 
Xc- Realistic cosmology obtains for E5i ~ 1/138. 

The cosmological dynamics of CQ is described in [14, 10]. For given E and 
Js = 0.05 we solve the flow equation (9). The "initial value" <5j = 5(x = rn) is tuned 
such that at present equals the value extracted from observation. The function 
5(x) in the effective action (6) is now fixed and we can derive the field equations for 
the evolution of gravity and the cosmon in a homogeneous and isotropic universe 
[10]. For simplicity we neglect the small influence of radiation and matter on the 
cosmon field equation which results from dZ F /d\nx ^ 0. For a solution of the field 
equations we start for x(t — 0) — m with rather arbitrary values of x(0) and Q^,(0) . 
The late time behavior will not depend on this. 

In table 1 we present a few characteristic cosmological quantities, namely the 
present Hubble parameter (h), the present fraction of dark energy, ^\ the present 
equation of state, \ the age of the universe, t^°\ the fraction of dark energy 

(Is) 

at last scattering, fi^ , the location of the third peak in the CMB anisotropies 
in angular momentum space, I3, as well as the normalization of the spectrum of 
density fluctuations, cr 8 , divided by the value which would obtain for a cosmological 
constant with the same fi^. Crossover quintessence with cosmological parameters 
in this range is well compatible [20] with the recent precision measurements of the 
CMB- anisotropies by WMAP [21]. 



E 


h 




(0) 


t (o) / 10 i<y 




^3 


/ (A) 


5 


0.66 


0.70 


-0.93 


1.37 


0.019 


796 


0.70 


12 


0.71 


0.73 


-0.99 


1.35 


0.0083 


794 


0.85 



Table 1: Characteristic cosmological quantities for two CQ models. 

We recall that smaller h and larger f^ "* shift I3 to larger values. We conclude that 
for large enough E our model is consistent with present cosmological observations. 

In a grand unified theory the value of the electromagnetic fine structure constant 
depends on the gauge coupling at the unification scale, ax, the ratio between the 
weak scale and the unification scale, Mw/Mx, and particular particle mass ratios 
like m n /Mw or m e /Mw for the nucleons and electrons. Since the nucleon mass is 
determined by the running of the strong gauge coupling similar dependencies arise 
for the ratio m n /Mx- 

For our model we assume that the dominant effect is due to the variation of «i 
[10] and consider the approximation of time-independent mass ratios Mx/M, Mw/m,, 
Mw/m e etc. The dependence of the fine structure constant on redshift can then 
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be related directly to the /^-function for the grand unified gauge coupling (7). For 
Aa em (z) = a em (z) - a em (0) one finds 



- — -2-Ao: x (z). (12) 



The same approximation yields 7 

A(m n /M) _ it Aa em _ ^ ^ Aa em 
(m n /M) lla em a em a em 

Grand unification links the time variation of a em with the time variation of m n /M 
or the effective Newton constant, G/G = 2d\n(m n / M) / dt = 78.2d em / a em . 

In terms of the variable x = In a = — ln(l + z) the time dependence of ax(x) 
obeys (cf. [10]) 

-dx~ = ^-frT = V3 ^ (1 + Wh) - (14) 

One observes that the z-dependence of Aa em is particularly weak in the region 
of large 5 and Wh close to — 1. This will be our explanation why R is consid- 
erably smaller than expected from a simple extrapolation linear in z. Crossover 
quintessence is precisely characterized by large 5 and Wh ~ — 1 in the region of small 
z, whereas y/1 + Wh/ VS is considerably larger for intermediate z ~ 1.5 — 3. One also 
observes that for constant +w h )/5) 1 ^ 2 and (3 a /a the dependence of Aa em (z) 

is logarithmic in 1 + z. This is crucial for large z as for nucleosynthesis, where Aa em 
turns always out to be much smaller than expected from a linear extrapolation in 
redshift or time. 

In eq. (14) the factor 1 + Wh is generic for all scalar models since it reflects the 
kinetic energy and therefore the rate of change of the scalar field. The substantial 
change in this factor is the main effect for a possible explanation of a small value of R 
by CQ. The factor (£\/<5) 1//2 is constant for small i\ in our model. It may decrease 
at late time if S grows large since i\ < 1 remains bounded. Finally, the details 
depend on the form of j3 a . In our model the fine structure constant would decrease 
with time for b 2 > and 6 6 = 0, at a rate that would not be observable for the value 
b 2 = 0.2 considered here. However, the influence of the change of 5 on the location 
of the fixed point ax,* leads to an increase of a em . As a consequence, the overall 
normalization of the time variation strongly depends on fe 6 . As a benchmark we use 
eq. (1). Therefore, for given b 2 and J a we fix such that Aa em (z = 2) = —7 - 10~ 6 . 
The shape of Aa em (z) depends now on b 2 and J a as well as on on E. In table 2 we 
show the values of Aa em at z = 0.13,0.45, 1100 and 10 10 for b 2 = 0.2, J a = 6. An 
investigation of other choices of b 2 and J a (in a range b 2 > 0.05) typically shows (for 
given E) a spread in the value of R by a factor rs 2. 

7 This relation differs from [22] where Mw/Mx is kept fixed instead of fixed m n /Mw in the 
present work. The quantitative difference is not very important, however. 
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E 


z = 0.13 


z = 0.45 


: = 1100 


z = 10 1U 


R 


V 


12 


-6 • 10~ 8 


-2.6 • 10" 7 


-6.5 • 10" 5 


-7.8- 10" 5 


0.0085 


5.1 • 10~ 14 


5 


-1.7 • 10" 7 


-8.2 • 10" 7 


-4.2 • 10" 5 


-5.4- 10" 5 


0.024 


4.4 • 10" 14 



Table 2: Time variation of the fine structure constant Aa em (z)/a em for various 
redshifts z. The last two columns show R (eq. (4)) and the differential acceleration 
V- 

Consistency with the Oklo bound (2) and the Re-decay bound (3) is found for 
large E. By looking at smaller values of E it also becomes clear that a precise 
knowledge of Aa em (z) could be used as a probe for the dynamics of quintessence. 
For example, it seems very difficult to achieve R < 1/50 for small values of E, quite 
independently of the precise form of (5 a . The present investigation confirms that a 
small value of R requires <C 1. This property of the equation of state seems 

to be favored by the cosmological tests (cf. table 1). We conclude that the bound 
(5) can be justified within a very broad class of models if one requires \R\ < 0.02. 
Without a substantial change in 1 + Wh between z = 2 and z = 0.13 it is very 
difficult to obtain a large change in da em /d\n(l + z) (cf. eq. (14)) as required for 
small \R\. Only for the large value E — 12 the functional dependence of Aa em (z) 
deviates substantially from a linear behavior in z in the range 0.5 < z < 2, with a 
strong "jump" of Aa em by a factor of 5 between z — 1 and z — 2. 

Let us next turn to the high values z ~ 10 10 characteristic for nucleosynthesis. 
In our scenario the dominant effect is a change in m n /M which modifies the clock 
[1], i.e. the rate of decrease of temperature in time units ~ m^ 1 - A decrease of 
m n relative to M is equivalent to an increase of M at fixed m n (and thereby to a 
decrease of Newton's constant). In turn, this decreases the Hubble parameter for 
a fixed temperature T (measured in units of m n ). Assuming that the primordial 
helium abundance Y p is in agreement with the value for time independent couplings 
within an uncertainty \AY P /Y P \ < 8- 10~ 3 and using A\nY p = |Aln(m n /M), yields 
the bounds |Aln(m n /M)| < 0.025 or 

,Aa em (z = 10 10 ) . „ , . 

| — -| < 0.64 ■ 10" 3 . (15) 

The effect of the variation of m n j M is much bigger than the influence of the variation 
of the proton-neutron mass difference due to electromagnetic effects. This explains 
why the bound (15) is more restrictive than several previous bounds [8], [9]. We infer 
from table 2 that our model meets the bound (15) as well as weaker bounds from 
the CMB [9]. It is striking that the values of Aa em /a em at nucleosynthesis are very 
far from an extrapolation with d In a em /dt = const, or even from d In a em /d In a = 
const.. 

In our approximation where the dominant field dependence of the various cou- 
plings arises from f3 a we can also estimate the size of the differential acceleration n 
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between two test bodies with equal mass but different composition. One finds [10] 



, = -6.10->f^)V'Afl I = -1.75.10-f^«=)' ,„, Aflz m (16) 

where for typical experimental conditions ARz = AZ/(Z + N) ~ 0.1. Here f3 a /ax 
and 5 have to be evaluated at z — 0. Results for rj are also displayed in table 2. We 
conclude that our CQ-model is compatible with the present experimental bounds [11] 
\r)\ < 3- 10~ 13 . We note that for a given value of d In a em / dz at z = (e.g. as inferred 
from Oklo) the differential acceleration is strongly enhanced for 1 + < 1. In 
addition, the effects of variations of mass ratios like M w /m n or m p /m e (that we have 
neglected here) can substantially enhance r\ [10]. Already at the present stage the 
tests of the equivalence principle can help to discriminate between different models! 
Further improvements of the accuracy of tests of the equivalence principle should 
either confirm or reject the interpretation of the QSO-results within our CQ-model. 

In summary, we have found models of crossover quintessence that can explain the 
QSO value for the time dependence of the fine structure constant and are neverthe- 
less compatible with all observational bounds on the time variation of couplings and 
tests of the equivalence principle. This explanation does not involve any particular 
cancellation of effects from the variation of different couplings. The same models do 
also very well with all present cosmological observations! 

What will be the tests of such models? First of all, the QSO result should 
be confirmed by new independent data. Hopefully this will allow for more precise 
restrictions on the shape of the function A \na em (z) which can be directly compared 
with the model predictions. Second, an improvement of the accuracy of the tests of 
the equivalence principle by an order of magnitude should lead to a direct detection 
of the new interaction mediated by the cosmon field. In view of its important role 
in present cosmology we may call this fifth force (besides gravity, electromagnetism, 
weak and strong interactions) "cosmo-dynamics" . Third, cosmological tests should 
distinguish crossover quintessence from a cosmological constant. 

The main signature in this respect is perhaps not the equation of state of dark 
energy in the most recent epoch (say for < z < 0.5). The value of in table 1 
shows that this quantity may be quite near to the value for a cosmological constant 
(i.e. = —1). Tests of Wh(z) at z > 1 by supernovae or gravitational lensing 

may be more sensitive. Indeed, if the QSO-results on the time variation of the fine 
structure constant are to be explained by quintessence this requires a substantial 
evolution of the cosmon field in the epoch near z = 2 (as compared to the rate of 
evolution today). We therefore expect that 1 + Wh(z) increases with increasing z. 
In turn, this implies that Vth(z) is substantially larger during structure formation 
and last scattering as compared to the value expected for a cosmological constant. 
Observations of structure formation and the CMB may be able to measure the 
value of a 8 /o"g A ' > and fij^ (table 1) with sufficient precision for a distinction from 
the cosmological constant. Finally, an improvement of the accuracy of atomic clocks 
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by three orders of magnitude could replace the bound (2) by an equivalent bound 
on d In a em / dz at z = within an experimentally well controlled setting. 



References 

[1] C. Wetterich, Nucl. Phys. B302 (1988) 668; Nucl. Phys. B302 (1988) 645; 
Astron. Astrophys. 30 (1995) 321 

[2] G. Dvali, M. Zaldarriaga, hep-ph/0108217; T. Chiba, K. Kohri, hep-ph/0111086 

[3] J. K. Webb et al., Phys. Rev. Lett. 87 (2001) 091301 

[4] J. -P. Uzan, hep-ph/0205340 

[5] A. I. Shylakhter, Nature 264 (1976) 340; 

T. Damour, F. Dyson, Nucl. Phys. B480 (1996) 37; 
Y. Fujii et al, Nucl. Phys.. B 573 (2000) 377 

[6] K. Olive et al., Phys. Rev. D66 (2002) 046007 

[7] C. Wetterich, hep-th/0210156 

[8] E. Kolb, M. Perry, T. Walker, Phys. Rev. D33 (1986) 869; 

L. Bergstrom, S. Iguri, H. Rubinstein, Phys. Rev. D60 (1999) 045005; 
K. Nollet, R. Lopez, Phys. Rev. D66 (2002) 063507; 
K. Ichikawa, M. Kawasaki, hep-ph/0203006 

[9] P. Avelino et al, Phys. Rev. D64 (2001) 103505; 
J. M. Martins et al, Phys. Rev. D66 (2002) 023505 

[10] C. Wetterich, hep-ph/0203266 

[11] S. Baessler et al, Phys. Rev. Lett. 83 (1999) 3585 

[12] B. Ratra, J. Peebles, Phys. Rev. D37 (1988) 3406 

[13] E. Copeland, A. Liddle, D. Wands, Phys. Rev. D57 (1998) 4686; 

R. Caldwell, R. Dave, P. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582 

[14] A. Hebecker, C. Wetterich, Phys. Lett. B497 (2001) 281 

[15] M. Birkel, S. Sarkar, Astropart. Phys. 6 (1997) 197; 
R. Bean, S. Hansen, A. Melchiorri, astro-ph/0104162 

[16] M. Doran, M. Lilley, J. Schwindt, C. Wetterich, ApJ. 559 (2001) 501 

[17] P. Ferreira, M. Joyce, Phys. Rev. D58 (1998) 1582 



9 



[18] M. Doran, J. Schwindt, C. Wetterich, Phys. Rev. D64 (2001) 123520 

[19] T. Damour, A. Polyakov, Nucl. Phys. B423 (1994) 532; 

T. Damour, F. Piazza, G. Veneziano, Phys. Rev. Lett. 89 (2002) 081601; Phys. 
Rev. D66 (2002) 046007 

[20] R. Caldwell, M. Doran, C. Mueller, G. Schaefer, C. Wetterich, astro-ph/0302505 

[21] C. L. Bennet et al., astro-ph/0302207; 
D. N. Spergel et al, astro-ph/0302209 

[22] X. Calmet, H. Fritzsch, Eur. Phys. J. C24 (2002) 639; Phys. Lett. B540 (2002) 
173; hep-ph/0211421; 

P. Langacker, G. Segre, M. Strassler, Phys. Lett. B528 (2002) 121 



10 



